Quantum Hall to charge-density- wave phase transitions in ABC-trilayer graphene 
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ABC-stacked trilayer graphene's chiral band structure results in three (n = 0, 1,2) Landau level 
orbitals with zero kinetic energy. This unique feature has important consequences on the interaction 
driven states of the 12-fold degenerate (including spin and valley) N = Landau level. In particular, 
at many filling factors vt = ±5, ±4, ±2, ±1 a quantum phase transition from a quantum Hall liquid 
state to a triangular charge density wave occurs as a function of the single-particle induced LL 
orbital splitting All- Experimental signatures of this phase transition are also discussed. 
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Introduction- When a group of partially filled Landau 
levels (LLs) with internal flavor indices are nearly degen- 
erate, electron-electron interactions can lead to sponta- 
neous broken symmetries which induce quasiparticle gaps 
and hence interaction driven integer quantum Hall (IQH) 
effects [lj. In semiconducting two-dimensional electron 
systems (2DES) where Zeeman splitting is negligible and 
the LLs consist of narrowly spaced doublets, interactions 
lead to spontaneous electron spin polarization Q accom- 
panied by a large exchange enhanced spin gap. More 
recently interaction driven broken symmetry states at in- 
teger fillings have been observed in both graphene d, Q 
and bilayer graphene 0, Q which in most cases exhibits 
a set of four-fold degenerate LLs. At many filling factors 
this four-fold degeneracy, which is due to the spin and 
valley degree of freedom, is spontaneously broken result- 
ing in charge gaps and hence incompressible quantum 
Hall (QH) states. Bilayer graphene's chiral bands lead 
to an additional degeneracy doubling of the zero-energy 
LL associated with the n = and n = 1 LL orbitals [7|, 
leading to an eight-fold degenerate LL. This additional 
LL degeneracy also results in interaction driven QH ef- 
fects due to an exchange induced gap between the LL 
pseudospins [1, Hj]. 

At relevant energy scales few layer graphene with dif- 
ferent stacking sequences, referred to as chiral 2DES 
(C2DES), exhibit electronic properties distinct from 
semiconducting 2DES [9( • The difference in the LL struc- 
ture of C2DES is most evident in the zero-energy LL: for 
certain stacking sequences (such as AB, ABC, ABCB, 
etc.) multiple LL orbitals are degenerate [l(J. The role 
of interactions on the ground state when multiple LL 
orbital degrees of freedom are degenerate remains rel- 
atively unexplored. In this Letter we show that, un- 
like bilayer graphene, at certain filling factors interac- 
tion induced QH states are unstable in ABC-stacked tri- 
layer graphene's 12-fold degenerate zero-energy LL. This 
is due to the presence of a degenerate set of triplet 
(n = 0,1,2) LL orbitals. While the QH liquid state is 
stable for v = —3, 3, 6 the nature of the ground states 
at v = ±5,±4,±2,±1 depends on the ratio of a single 
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FIG. 1. (Color online) Schematic representation of the inter- 
action induced energy gaps as a function the magnetic field 
for interaction driven states in ABC-trilayer graphene's 12- 
fold degenerate LL, the numbers represent the Hall conduc- 
tivity in units of e 2 /h. For (a) All < A L C L the presence 
of insulating triangular CDW states for v — ±5, ±4, ±2, ±1 
results in Hall plateaus only at a xy = ±6, ±3, 0(e 2 /h). For 
(b) All > A L C L interaction driven QH states appear at all 
intermediate integer fillings with the sequence determined by 
Hund's rule behavior (see text for details). 



particle induced LL orbital splitting All and electron- 
electron interactions e 2 /els (Jb — y/hcfeB is the mag- 
netic length). In particular, a quantum phase transition 
occurs from a triangular charge density wave (CDW) to 
a translationally invariant QH liquid state as All is in- 
creased beyond a critical value A^[l2j]. This transi- 
tion is first order in the liquid-solid universality class. 
The crystal state exhibits coherence between LL orbitals 
which leads to vortex textures of in-plane electric dipoles 
at each lattice site. For moderately disordered samples 
this phase transition should be characterized by a Hall 
plateau transition as a function of All at a fixed filling 
factor v. Another consequence is that the Hall conductiv- 
ity for All < ^ll wm correspond to the adjacent inter- 
action driven IQH plateau, for example, at v = —5, —4, 
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a xy = — 6e 2 /h (as shown in Fig.Q^). 

ABC-trilayer graphene- ABC-stacked graphene trilay- 
ers are chiral generalization of monolayer and Bernal- 
stacked bilayer graphene. ABC-stacked trilayer is a par- 
ticular stacking of three graphene layers each with in- 
equivalent sublattices A4 and Bi arranged in following 
sequence: one of the two-carbon atom sites in both 
the top and bottom layer Bi(A-z) has a different near- 
neighbor carbon atom site in the middle layer ^2(^2), 
which leaves one-carbon atom site in the top and bot- 
tom layers Ai(B^) without a near- neighbor in the middle 
layer. Interlayer hopping on adjacent layer near-neighbor 
carbon atom sites leads to the formation of high-energy 
dimer bands which push the electron energy away form 
the Fermi surface leaving one low-energy sublattice site 
per 7r-carbon orbital in the outermost layers. Because 
the hopping between the low-energy sites via intermedi- 
ate high-energy states is a three-step process, it turns out 
that in the presence of an external magnetic field three 
LL orbitals (n = 0,1,2) have zero kinetic energy. The 
wavefunctions of the opposite K(K') valleys in the zero- 
energy LL are localized in the top(bottom) layers, hence 
layer and valley labels are equivalent. This simplified 
two-band model for the ABC-stacked graphene trilayer 
which includes only nearest-neighbor interlayer and in- 
tralayer hopping terms is valid at large external magnetic 
fields Iff . The effects of an external potential difference 
between the top and bottom layers Ay along with remote 
next-nearest interlayer hopping terms on the zero-energy 
LL, not included in our simplified model, can be incorpo- 
rated by introducing a symmetric LL gap All [3 be- 
tween the degenerate zero-energy triplet set in = 0, 1, 2) 
of LL orbitals. The effect of these terms on All will be 
reported elsewhere [la ]. 

Mean field theory- Trilayer graphene's zero-energy LL 
is a direct product of two S — 1/2 doublets corresponding 
to spin and valley (layer) pseudospin along with S = 1 
LL orbital pseudospin which is responsible for the new 
physics discussed in this Letter. Working in the Landau 
gauge A = (0, Bx, 0), the Hartree-Fock (HF) mean-field 
Hamiltonian, projected onto the zero-energy LL of ABC- 
stacked trilayer, consists of single particle pseudospin 
splitting terms along with direct and exchange interac- 
tion contributions: Hhf = H sp + Hint, with 
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(repeated indices are summed over) where n., = 0, 1, 2 are 
the LL orbital indices, a, (3 — t(b) label the top (bottom) 
layers (K(K') valleys) and r,s = 1(-1) label the t (i) 
spins, = A/(2tt1 2 b ) is the degeneracy of a single LL 



and momentum space guiding center density is , 

KtpM) = E C ix aS Cn'X' re-^ X ' + ^S qylltX _ x , 

x,x> 

(2) 

where c\ lXas {c n xas) are the respective n th LL creation 
and annihilation operators at a given guiding center X 
in the Landau gauge with layer and spin index a, s. The 
effects of the single-particle pseudospin splitting field are 
encoded in 
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where p n as = A%£ as (q = 0) = Ex( c lx as c nX as ) is the 
total guiding center density for the n th LL with valley 
and spin index a, s and All is the LL splitting induced 
by the remote interlayer hopping terms and the external 
potential difference Ay with £ t (M = 1(— 1), Az = g^sB 
is the Zeeman splitting and all energies are measured in 
units of e 2 /elB- In CQ) the Hartree-field captures the 
valley (layer) dependent electrostatic contributions: 
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and the exchange contribution is captured by: 
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[la] is 



the Fourier transform of electron-electron interactions in 
the same and different layers (valleys) and d = 0.667nm 
is the interlayer separation. The form factors F nn i (q) 
capture the character of the three different LL orbitals 
and are given by: 



Fun' (q) = \/ — r 
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for n' < n, where L^(x) is the generalized Laguerre poly- 
nomial, 9 q — ta,n~ 1 (q y /q x ) and the form factors satisfy 

Fnn'(q) = [^n'nC-q)]*- 

First we look for translationally invariant solutions of 
the density matrix (A^™ a , s , (q = 0)). The density matrix 
must be determined by self-consistently occupying the 
lowest energy eigenvectors of Hhf- For a uniform den- 
sity order parameter the only contribution to the Hartree 
term is Eh — (e 2 /e^)(d/2/s), which captures the charg- 
ing capacitance energy of the outermost layers of ABC- 
trilayers. The filling order of the 12-fold degenerate zero- 
energy LL proceeds in integer increments starting from 
the filling factor v = — 6 and follows a Hund's rule be- 
havior: first maximize spin polarization, then maximize 
the layer pseudospin polarization to the greatest possible 
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FIG. 2. (Color online) Collective mode dispersion uj q for 
v = —5,-2,1,4 in units of interaction strength e 2 /e/s = 
11.2(B[T]) 1/2 meV as a function of qls at different values of 
the LL splitting All at a magnetic field of B = 10 Tesla. 
Charge density wave instability is accompanied by a soft- 
ening of the magneto-roton minima for All < ^l*l- The 
gapped mode corresponds to the hard-direction in LL pseu- 
dospin space and is not essential for our discussion (see text 
for details). 



extent and finally maximize the LL polarization to the 
greatest extent allowed by the first two rules (see Fig.[Tj3). 

For balanced ABC-trilayers (Ay = 0) at v = —3, 3 
the valley(layer) pseudospin exhibits spontaneous inter- 
valley(layer) coherent state similar to that observed in 
conventional 2DEGs [Hj]. The v — state is spin po- 
larized as the exchange gap associated with spin polar- 
ization is larger than the exchange gap associated wiht 
the inter-layer coherent state. However, since the en- 
ergy scales associated with Zeeman splitting and the ratio 
cL/Ib are small compared to typical interaction energies, 
the valley and spin exchange gaps compete closely which 
could alter the Hund's rule behavior (for example leading 
to density wave states rather than spin-polarized states 
at v = 0). Another consequence of a small value of (1/Ib 
is that for a small value of A v (> (vr/8) 1/2 (d// B ) 2 e 2 /£/ B ) 
the inter-layer coherent state switches to a valley polar- 
ized state. 

Hund's rule dictate that the topmost LLs at v = 
—5,-2,1,4 have n = LL orbital occupation whereas 
at v = —4,-1,2,5 the n = 1 LL orbital is occupied. 
This preference of LL orbital occupation is associated 
with a better exchange energy of the more localized 
n = LL orbital wavefunction along with the absence 
of nodes The new physics discussed in this Letter 
is primarily associated with the LL pseudospin ordering 
at ^7^—3,0,3,6. In the next section we show that for 
All < ^ ne translationally invariant uniform liq- 

uid QH states at v = —5,-2,1,4 are unstable to the 
LL orbital fluctuations. This instability is due to strong 



charge fluctuations which lead to a softening of the col- 
lective mode dispersion as a function of All- 

Charge density wave instability — We now focus on the 
collective excitations of a filled QH liquid at filling fac- 
tors v = —5,-2,1,4, freezing spin, and valley (layer) 
but keeping the (n = 0, 1, 2) LL orbital degrees of free- 
dom. To further simplify the discussion and concen- 
trate on the LL pseudospin fluctuations, we also assume 
SU(4)-invariant spin and valley (layer) independent inter- 
actions. [l9| In this case the SU(4)-flavor index is sponta- 
neously broken and quantized LL pseudospin fluctuation 
consists of pseudospinor rotations which mix n = with 
n = 1 and n = 2 LL orbitals leading to two branches of 
collective mode excitations. The collective mode excita- 
tions which we evaluate below are primarily influenced 
by the orbital dependence of the microscopic Hamilto- 
nian fl} . A detailed account of the collective excitations 
for d 7^ will be presented elsewhere [l5[ ■ 

One physically transparent way of performing these 
collective mode calculations is to construct a fluctuation 
action in which each transition has canonically conjugate 
density p and phase <p components corresponding to the 
real and imaginary parts of the final state component in 
the fluctuating spinor. The fluctuation action 

S\P,v]=Sb\p,v\- J du J d 2 q£[p,tp], (7) 

contains a Berry phase term Sb [2Cj| 

S B [P,(P\= J du} J d 2 q ipl q XV q - ^ q 2? f pq J (8) 

where V = -iuhx2 and p q = (pi, q , P2,q) with p n ^ cor- 
responding to the charge fluctuation of the n = to the 
n th orbital, with a similar definition for the phase vari- 
able (fq. The energy functional, £ [p, <p] is closely related 
to the discussion in the preceding section 
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q and Aq are 2x2 Hcrmitian matrices which cap- 
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ture the energy cost of small pseudospinor fluctuations 
and can be evaluated explicitly. The matrix elements of 
T q can A 9 can be expressed in terms of direct and ex- 
change matrix elements: 



Hi = 

Kj = 



(HSU + (-l) i+i i%o(|q|,0) - X mj (q, 0), 
(-l)^^ OjO (|q|,0)-X OOiJ (q,0), 



(10) 



where a x = A LL + l/2^/ir~/2 and a 2 = 2A LL + 5/8\/7r72 
are wavevectors independent contributions measured in 
units of e 2 /eIb- 

The collective mode dispersion plotted in Fig. [2] indi- 
cates that the uniform QH liquid state is unstable for 
A LL < A ( * } = 0.015e 2 /(d B ). The gapped mode in 
Fig. [5] which corresponds to some wavevector dependent 
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linear combination of the orbital pseudospin fluctuations 
(Pn.qj Vn.q) reflects the difference in the exchange energy 
cost of the LL orbitals. The "soft-mode" indicates the in- 
stability of the uniform QH liquid at qol B ~ 2.37 which 
is due to the charge fluctuations of the QH liquid as we 
show next. 

In order to analyze this instability and its associated 
consequences it is convenient to expand the energy func- 
tional of the "soft-mode" density p_, q near q ~ q , 

£ [p_, q ] = [a(A LL ) + a&LL)(q - qo) 2 ] p 2 _, q (H) 

+ /?(A LL )p_ qi p_ q2 p_ q3< 5(qi + q 2 , q 3 ) H 

where a{A LL ) = 1.24(A LL - A L * ] L ) with £(A LL ) = 
0.29 - 0.14A LL and /3{A LL ) is the coefficient of the 
third order term which is not prohibited by any symme- 
tries. For /3 = (ITT1) exhibits a Brazovskii-like instability 
leading to a fluctuation-induced isotropic-smectic transi- 
tion to a unidirectional LL-pseudospin density wave pat- 
tern [21] . However, the presence of a finite j3 preempts 
this transition to a triangular two-dimensional CDW (or 
crystal) [22j with LL orbital pseudospin textures. 

This instability should be present whenever (n > 
2) LL orbitals are degenerate. The effective inter- 
action within the projected n th LL orbital v n (q) = 
v q [L n (q 2 lg/2)] 2 e~ q ~ l B/ 2 has minimas corresponding to 
zeros of the Laguerre polynomials. Due to this the 
Hartree contribution vanishes at these finite-g wavevec- 
tors and since the exchange energy is positive for all 
wavevectors, E B f becomes negative 23]. The near de- 
generacy of n = 0, 1, 2 LL orbitals allows electrons access 
to lower energy by forming coherent superposition of LL 
orbitals for a range of finite-g wavevectors, leading to a 
CDW instability. As All is increased from below this 
gain in energy competes with the single particle contri- 
bution and the uniform QH state becomes energetically 
more favorable resulting in a transition. 

Triangular Charge Density Wave (All < A LL ) — 
For All < A [ [ ] L with A [ [ ] L = 0.095e 2 /el B > A [ * ] L at 
v = —5,-2,1,4, the HF solution with the lowest en- 
ergy is a triangular crystal of electrons with one electron 
per lattice site as shown in Fig. |3] |24j. In this crystal 
state, all three orbitals n = 0, 1, 2 are occupied in contrast 
with the QH liquid state where only n — is fully occu- 
pied. As All ~> A L c) L fr om below, the total population 
(Ai,i) + (A 2 ,2) of the higher-energy orbitals decreases but 
stays finite until A^ where it drops discontinuously to 
zero indicating a first order transition to the liquid state 
with (Ao,o) = 1. At each crystal site, the electron is in a 
superposition of the three orbital states n = 0, 1, 2. Since 
the wave function profile of the three orbitals is different, 
this superposition leads to the triangular CDW pattern 
in the total electron density n (r) shown in Fig. [3] (a). 
When two orbital are filled (at v — —4, — 1, 2, 5 for exam- 
ple), the ground state at small All is a crystal of holes 
instead of electrons (see Fi^3] (b) ). 
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FIG. 3. (Color online) Electronic density and orientation of 
the in-plane electric dipole field for the triangular crystal at 
filling factor (a) v = -5, -2, 1, 4 and (b) v = -4, -1, 2, 5. The 
electrons for (a) and holes for (b) crystallize in a triangular 

pattern with a lattice constant a = <J 4tv/\/$Ib- The crystal 
state is characterized by a superposition of three LL orbitals 
n — 0, 1, 2 pseudospin which manifests itself as a periodic 
density modulation n (r) (see text for details). 

The CDW state is characterized by a periodic mod- 
ulation of (A„ iin (q)) (n 7^ m) defined in Eq. ([2]). This 
leads to the presence of an in-plane electric dipole field (a 
similar situation also occurs in bilayer graphene[25j]): 

dr =x , y (q) = -V2d B e-i 2l */ A (q)) (12) 

-2d B e-W (1 - g 2 /|/4) ( P P (q)), 

where pi"' m) = (A„. m + A rn , n ) /2, p { y ,m) = 
(A n m — A m n ) /2i (with n > to), the energy of 
the CDW state in the presence of an external electric 
field is - / dr (d (r) • E (r)) . Fig. EUa)-(b) also show the 
orientation of the in-plane electric dipole field around 
each crystal site which exhibits vortex textures in 
which the dipole vectors rotate by 2ir. The topological 
properties of this crystal state will be characterized 
elsewhere (l5| . 

The degeneracy of n = 2 LL orbital is crucial for the 
existence of the crystal state. This situation is somewhat 
reminiscent of the various crystal phases that appear 
in conventional 2DES at fractional fillings of higher 
LLs n > 2 [2^, [2(|. The surprising result here is that 
similar electronic crystal phases appear at total integer 
filling factors in ABC-stacked trilayer graphene and are 
sensitive to the competition of interactions and single 
particle terms induced by remote hopping and layer 
imbalance. This triangular CDW state would likely 
be pinned by disorder leading to insulating behavior 
of the uppermost LL thus exhibiting quantized Hall 
conductivity at a value corresponding to the adjacent 
interaction driven IQH plateau (see Fig. QJi). The 
triangular CDW state should also be visible through its 
pinning mode in microwave absorption experiments. 
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